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1 INTRODUCTION 



The recently catagorised L-dwarfs (Kirkpatrick et al. 199£) 
have rapidly become an important part of stellar astro- 
physics. L-dwarfs lie below M-dwarfs in terms of their tem- 
peratures, and their spectra are dominated by molecules 
such as TiO, FeH and CrH. Dust formation also plays a 
very significant role in the atmospheric physics of early L- 
dwarfs (by L5-7, the major dust species have begun to settle 
below the photosphere). The physics of dust formation is a 
complicated subject, but one which is critical to our under- 
standing of L dwarfs. Variability studies give us a tool to 
probe the sub-global differences in the composition of the 
atmosphere. Other processes, such as the flaring seen in M 
dwarfs, may also contribute to variability in L dwarfs. 

In this paper we report time-series photometry of four 
field L dwarfs. Two objects, 2MASS1108307+683017 and 
2MASS0746425+200032 are found to be variable at the level 
of 0.012 mag (5a) and 0.007 mag (6.5cr) respectively on the 
timescale of hours. We did not detect any variability above 
the noise level of ~0.02mag (la) from the other two targets 
(2MASS 1146345+223053 and DENIS0909571-065806). 



2 OBSERVATIONS AND REDUCTION 

Data were obtained with PFCam on the 3-m Shane tele- 
scope at Lick Observatory on the nights of 2002 January 
21 and 23 UT. Conditions were photometric on the first 
night, but patchy cloud interupted observations on the 23rd. 
PFCam has a 2048x2048 CCD array with a pixel scale of 
0.296" /pixel. The full field of view, 10.1' x 10.1', is more than 
sufficent for our purpose (high cadence, high precision differ- 
ential photometry). We therefore used only a small section 
of the chip to improve read-out times. Chip windows were 
choosen to select the target, and 5 or 6 comparison stars of 
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similar magnitude. The chip was also binned where possible 
to further reduce overheads, such that we collected images 
at 2 minute intervals on average. All observations were made 
through the I-band filter where the photon count from the 
extremely red L dwarfs is maximised. Calibration frames, 
consisting of bias frames, dome and sky flats were obtained 
for all chip configurations. A log of observations, including 
window sizes and chip binning, is given in table ^. 

Data were reduced with IRAFQ Ten bias frames were 
median combined, and this was then subtracted from all 
data and flat field frames. An average normalised flat field 
was then constructed from sky flats. We divided the science 
data by the appropriate flat field. This provided images flat- 
tened to the 0.25% level. Fringing is visible at the ~1% level 
on the scale of ~10". We could not remove this, as we did 
not obtain enough pointings to build a fringe map. 

During the first night's observing, we discovered PF- 
Cam's shutter was working incorrectly, with only one blade 
of the traveling shutter system functioning. As the travel 
time of the blade is on the order of the milliseconds, a part 
m 10 6 for most science frames, this contributes neglibly to 
the photometric errors. The effect is however apparent in 
short (0.5-4s) dome flats. For these frames, we fitted the 
gradient with a linear function, and subtracted the fit. 

In some science frames, the peak counts in a few pixels 
approached the full- well depth of the CCD (65536 ADU). We 
obtained data to test the linearity of the chip, and found it 
to be linear to ~60000 ADU. Any frames with peak counts 
above this level were removed from the analysis below. 

Photometry was performed with the APPHOT package 
within IRAF. We determined the FWHM for each frame 



1 The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by NOAO, which is operated by AURA, Inc., under con- 
tract to the NSF. 



2 F.J. Clarke, B.R. Oppenheimer and C.G. Tinney 



within an observing sequence, and picked an aperture ~1.5x 
larger than the average FWHM for that sequence. This aper- 
ture was used on all frames within that sequence. Magni- 
tudes and errors were calculated from t he fluxes given by 



APPHOT in the manner described by Everett & Howell 



(2001)[]including an estimate of informal errors from fring- 
ing and imperfect flatfield correction. An ensemble of bright 
stars were combined to give a zero-point magnitude for each 
frame (m ons ), which was then subtracted from the mea- 
sured magnitude of the target (m tar ) to build a differential 
lightcurve. Hence, an increase in 8m (m tar -m cns ) represents a 
dimming of the object. A selection of comparison stars (not 
included in the ensemble stars) were also analysed to ensure 
any variability detected is intrinsic to the target. In no case 
did any of the ensemble stars display any variability. 



3 RESULTS 



3.1 2MASSJ1108307+683017 



2M1108+6830 is an LI dwarf. |Gizis et al. (2000)| report an 
H-q equivalent width of 7.8A, which they suggest indicates 
it is possibly an older member of the stellar population, 
rather than a youthful brown dwarf. We observed 2M1108 
for 4.5 hours with a sampling rate of 2 minutes (100s in- 
tegrations). To obtain sufhcent ensemble stars (at least 3), 
we used a field of 5.l'x5.l' with 2x2 pixel binning (Fig- 
ure Q). Figure ^ shows the resulting I band lightcurves for 
2M1108 and a comparison star (which was not part of en- 
semble stars). It is clear from this figure that 2M1108 be- 
comes ~0. 012 mag fainter during the last 1.5 hours of obser- 
vation, whilst the comparison star stays constant to within 
±0.003 mag throughout the observation. The standard devi- 
ation of the binned comparison lightcurve is 0.0025 mag. Our 
detection of variability in the binned lightcurve is therefore 
significant at the 5a level. There may be a slight brighten- 
ing in 2M1108 just before it fades (AJD0=6.O3), but this is 
detected only at the 2a level due to an increase in noise in 
this region. 

Although the comparison star shows no variability, we 
must consider the possibility that 2M1 108's variability is 



caused by second order colour effects (see Young et al. (1991) 



for a full discussion). The L dwarfs have much redder spec- 
tral energy distributions than the ensemble stars, and hence 
have a longer effective wavelength in the I filter. This means 
the L dwarfs have a smaller effective extinction coefficent, 
which would manifest itself as an apparent brightening in 
the differential lightcurve with increased airmass. However, 
the data exhibit the opposite behaviour, since the observa- 
tion started at an airmass of 1.15 and ended at 1.3. The 
detected dimming of 2M1108 we see is therefore unlikely to 
be caused by airmass effects. 

The H-a emission from 2M1108 implies chromospheric 
activity. However, it is not clear if the variability we de- 
tect is related to this. The variability may alternatively be 
due to photospheric features (magnetic spots or dust clouds) 
which modulate the brightness as the object rotates. There 
is no periodic signal in our data, but our temporal sampling 



2 AJD = HJD-2452290.0 




HJD - 3452290.0 

Figure 2. I band differential lightcurve of 2M1108307+683017 
(upper panel) and the comparison star indicated in Figure |lj 
(lower panel). Individual measurements are shown with narrow 
lines, with the lightcurve averaged into 20 minute bins shown in 
thick crosses. 2M1108 shows a dimming (increasing 5m) of ap- 
proximately 0.012 mag over ~2 hrs. The comparison star stays 
constant to ±0.003 mag throughout the night. Both lightcurves 
have been mean subtracted. 



clearly does not cover a sufficent range t o exclude this pos - 
sibility. Based on wsini measurements (Basri et al. 2000), 
we expect the typical rotation period of an L dwarf to be on 
the order of 6-12 hours. 



3.2 2MASSJ0746425+200032 

2M0746+200 is infact a clo se (0.22") binary of roughly 
equal mass (Reid et al. 2001). We do not resolve the bi- 



nary in our observations, and therefore refer to the whole 
system as 2M0746 rather an either individual component. 

We discuss how this may effect our results later. 

2M0746 has a spectral type of L0.5. [Schweitzer et al 



(2001) detect no Li in its spectra, and therefore estimate a 
minimum mass of ~O.O6M0for each component. Due partly 
to its binary nature, 2M0746 is one of the brightest known 
L dwarf systems, having an I band magnitude of 15.2. It 
also lies in a rich stellar field, making it an ideal target for 
accurate differential photometry. We observed this target 
for nearly 6.5 hours with a sampling rate of 1-2 minutes. 
Changes in seeing required us to adapt our observing strat- 
egy (to avoid saturating the target). The majority of the us- 
able images measure 2.2' x 2.2', and contain the target and 
4 bright comparison stars (Figure ^). These stars are also 
visible in all previous configurations used on this object. 

The derived I band lightcurves are shown in Figure ^. 
2M0746 shows variability in two of the three configura- 
tions. In the second configuration (AJD=5.64-5.70), 2M0746 
brightens (<5m decreases) by ~0.01 mag while the compari- 
son stays constant to within 0.005 mag. However, due to the 
increased noise in this configuration, this is only signficant 
at the 2a level. During the third configuration (AJD=5.70- 
5.89), 2M0746 fades by ~0.007mag around AJD=5.78 be- 
fore returning to its original brightness. The comparison ob- 
ject is constant to ~0.002 mag, and its binned lightcurve has 
a standard deviation of 0.0011 mag. Our detection of vari- 
ability in the third configuration is therefore significant at 
the 6.5a level. 

It is possible that the brightening we see from 5.66 to 
5.7 is analogous to the brightening from 5.78 to 5.82. This 
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Object 1 


Date 


Start UT 


End UT 


# Frames 


Exp time 


CCD confif 


;uation 




2002 Jan 








w 


Window size 


binning 


2M0746425+200032 


21 


02:48 


03:33 


28 


45 


5.1'x5.1' 


2x2 


2M0746425+200032 


21 


03:33 


04:55 


51 


45-30 


3.7'x3.7' 


2x2 


2M0746425+200032 


21 


04:55 


09:22 


121 


120-60 


2.2'x2.2' 


lxl 


2M1108307+683017 


21 


09:52 


14:23 


119 


100 


5.1'x5.1' 


2x2 


D0909571-065806 


23 


07:51 


08:20 


6 


300 


2.2' X 2.2' 


2x2 


D0909571-065806 


23 


09:53 


11:57 


26 


300-200 


2.2'x2.2' 


2x2 


2M1146345+223053 


23 


12:11 


14:20 


22 


300 


2.2'x2.2' 


2x2 



Table 1. Log of observations o btained with PFCam o n the 3-m Shane Telescope at Lick Observatory. Ob ject names starting with 2M 
come from the 2MASS survey jKirkpatrick ct al. 1999 ). D0909-0658 was discovered by the DENIS survey ( |Delfosse et al. 1997 ). 
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Figure 1. Finder chart for all four targets, indicating the position of the target and the comparison stars (marked with arrows). Note 
the difference in scale of the 2M1 108+6830 field. 
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5.6b 5.7 5.75 5.8 

HJD - 2452290.0 



5.85 5.9 



Figure 3. I band lightcurve of 2M0746425+200032 (upper panel) 
and an apparently invariant comparison star (lower panel). In- 
dividual measurements are shown as narrow lines, with the 
lightcurve averaged into 20 minute bins shown with thick crosses. 
2M0746 clearly changes in brightness (dims) by ~0. 007 mag on 
the timescale of 3 hours while the comparison star stays constant 
to ~0. 002 mag. The dashed vertical lines at times ii.64 and 5.70 
indicate the changes in chip configuration (Table ph. Each con- 
figuration has a different zeropoint, and all three have been indi- 
vidually mean subtracted. Differences between the configurations 
are therefore not significant. 



may suggest that we are seeing rotational modulation of the 
lightcurve on a period of ^3 hours. If this were the case, 
we would expect to see the next minima at AJD~5.9 - just 
after the end of our observations. Three hours would be a 
fast, but not unprecedented rotational period for a brown 
dwarf. For example, Kelu-1 has a ro tational period of 1.8 
hours flClarke, Tinney fc Covey 2002| ). 

Unfortunately, the binary nature of 2M0746 complicates 
matters significantly. As we are observing an unresolved bi- 
nary, we are essentially sampling the combined lightcurve 
of two objects. It is not therefore possible to attribute spe- 
cific events in the lightcurve to either object. In addition, 
the binarity will tend to "wash-out" variability from either 
component of the binary; for example, 1% variability from 
one component will result in only 0.5% measured variability 
from the system. 

2M0746 is clearly an interesting target for future vari- 
ability observations, especially if the binary can be resolved. 



3.3 DENIS0909571-065806 



Martin et al. (1999) classify D0909-0658 as an L0 dwarf. 
Figure ^ identifies it along with the comparison star whose 
lightcurve is plotted in Figure ^j. Seeing on this night was ~2 
arcsec poorer than the first night, requiring longer (300s) ex- 
posures to obtain the necessary signal to noise ratio. D0909- 
0658 was observed in two batches, seperated by ~90 min- 
utes. Figure ^| shows no significant indication of variability 
either between or within the two sequences. We therefore 
place an upper limit on the variability of D0909-0658 of 
(Tj <0.02mag (1 sigma) on the timescale of 0.5-4 hours. 



3.4 2MASS1146345+223053 

2M1146+2230 is actually a close equal brigh t ness binary 
with a separ ation of 0.3" (Koerner et al. 199£). Schweitzer 
et al. (2001' have detected Li in the combined spectrum of 
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7.9 7.95 
HJD - 2452290.0 



Figure 4. Lightcurve of D0909571-065806 (upper panel) and the 
comparison star indicated in Figure |l] (lower panel). The is no 
significant variability from the target, with the scatter in both 
lightcurves at the level of 0^=0.02 mag. 
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HJD 



8.06 
2452290.0 



Figure 5. Lightcurve of 2M1146345+223053 (upper panel) and 
the comparison star (lower panel) denoted in figure [l| Both 
lightcurves are consistent with non-variability within the errors. 
This gives an upper limit of Oi <0.02mag for the variability of 
2M1146 over ~2 hours. 



2M1146, confirming its status as a brown dwarf binary with 
each component having M<0.06Mq. In addition, an earlier 
t ype background star i s located 1" away from the binary 
( Kirkpatrick et al. 1999 ) . The seeing during our observations 
(3 arcsec) did not allow us to resolve the binary or to detect 
the background star. 



Bailer-Jones & Mundt (2001) have previously made a 
marginal detection of variability in this object at the level 
of ~0.015 mag. Poor seeing did not allow us to reach this 
level of accuracy. The lightcurve (Figure ^ shows 2M1146 
is constistent with being non-variable at the noise level of 
~0.02mag (1 sigma) during our observations. 



4 CONCLUSIONS 

We have detected variability in time series I band ob- 
servations of two field L dwarfs (2M0746425+200032 and 
2M1108307+683017), and placed upper limits on the vari- 
ability of two more (2M1146345+223053 and D0909571- 
065806). The cause of variability is unclear, but it may be 
due to imhomogenous structures within the photosphere. In 
the case of 2M0746+2000, weak evidence exists for a peri- 
odicity of ~3 hrs. This could be explained by a photospheric 
feature and a 3 hour rotation period. 



L dwarf variability 5 



In now seems clear that a significant fraction of, if not 
all, L-dwarfs are variable at the 1-2% level, and that very 



few have much larger variability ( 


Clarke, Tinncy & Covey 


2002, 


t 


ailer- Jones & Mundt 2001 


, Martin, Zapatero Oso- 


rio & I 


ichto 2001). In this study, both stars for which wc 



obtained photometry better than 1% exhibited variability. 
Future surveys for variability in L-dwarfs require a photo- 
metric precision better than 0.5%. 
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